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Abstract 

We develop a field theoretical approach based on the temporary basis description as a tool to 
investigate the transmission properties of a time-driven quantum device. It employs a perturbative 
scheme for the calculation of the transmission of a monochromatic beam of particles through the 
time-dependent set-up. The main advantage of the proposed treatment is that it permits the use 
of the particle picture for the calculation of the scattering matrix and the transmission coefficient. 
Therefore the elementary physical processes contributing to the transmission can be identified and 
interpreted in a transparent way. We apply the method to the simple but prototype problem of 
transmission through an one-dimensional oscillating delta potential and we demonstrate how it 
enables a deep understanding of the underlying physical processes. 
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I. INTRODUCTION 



The transmission of quantum particles through a time-dependent potential has been the 
subject of extensive studies in the last three decades (see [l] and references therein). The 
main goal in these studies is to classify and understand the mechanisms of quantum tun- 
nelling in such a potential. Typical examples of processes involving transmission through 
time dependent devices are the tunnelling of a test particle from a metastable state [2] us- 
ing instanton techniques, the scattering off a dissipative environment and the tunneling 



through a time- modulated barrier j^. Particularly in j4| it was shown that at low mod- 
ulation frequencies the traversing particle sees a static barrier and at high frequencies the 
particle tunnels through the time-averaged potential. It was also demonstrated that inelas- 
tic processes can occur, where the tunneling particle gains or loses energy quanta from the 
modulation field. In the same work a fundamental question was raised concerning the defi- 
nition of the traversal time through such a fluctuating device, an issue which is still under 
debate jsl. Recently novel perspectives for the study of driven quantum systems emerged. In 
particular it has been realized that there is the possibility to exploit phenomena occurring 
in these systems for the development of novel technology such as the design of quantum 
pumps, i.e. devices capable to create quantum directed transport through periodic external 



driving js], or the control of quantum devices for information processing (7|. 

One of the simplest examples for studying the transmission in a time-dependent envi- 
ronment is the tunnelling through a delta-barrier with a harmonically oscillating coupling 
which has been considered by several authors js iO]. Although the device is very simple, the 
presence of the time-dependent coupling leads to interesting phenomena, observable in the 
transmission coefficient, such as Fano resonances, threshold enhancements related to side- 
band modes and their interplay jol. More complicated setups, involving several oscillating 
delta barriers, have been investigated in the context of quantum pumps In this case 
dephasing effects may also occur, influencing the appearance of resonances jisj. In jl^ an in- 
finite periodic chain of delta-barriers with harmonically oscillating strength has been studied 
demonstrating the modification of the conductance zones through the time-dependence. 

If the driving of the potential barrier is periodic the standard approach used to solve the 



quantum dynamics is based on Floquet theory 



14|] which is very well suited for the straight- 



forward calculation of the transmission properties, discriminating between the contributions 



from elastic and inelastic channels. In the Floquet treatment the considered problem is 
first rendered time-independent and then solved using standard numerical techniques. This 
represents certainly an advantage from the computational point of view. However, it is 
conceptually a handicap since it does not allow to isolate contributions or identify mecha- 
nisms leading to a specific dynamical behaviour that allow for a better understanding of the 
underlying physics. 

Path integral methods, on the other hand side, have been commonly used for the study 
of the quantum dynamics in time-dependent potentials [15] but they show usually a slow 
convergence. Thus, despite of a few exceptional cases where an analytical solution can 
be obtained an efficient general purpose path integral treatment of time- dependent 
potentials is still lacking. Progress in this direction has recently been achieved by introduc- 
ing a rapidly converging scheme in the framework of a high-order short-time expansion of 



transition amplitudes in time- dependent potentials [17|. This approach holds for a general 
time-dependence of the potential. The disadvantage of this approach is that it can only 
be applied to smooth potentials. This is a relevant point in view of the fact that previous 
works have mostly been employing either the 5 or a rectangular barrier for the study of the 
time- dependent systems. 

In the present work we develop a field theoretical approach to the investigation of the 
quantum dynamics in time-dependent potentials. The method is suitable for studying any 
kind of potentials (also non-smooth) as well as external driving (periodic or not). In addi- 
tion it allows for a classification and decomposition of the quantum dynamics in terms of 
fundamental processes which represent the skeleton of the quantum evolution in these sys- 
tems. Particularly we reveal the role of virtual "multi-photon" processes which determine 
to a large extent details of the resonant structures in the dependence of the transmission 
coefficient on the incoming energy. Within our treatment these sub-processes can be clearly 
distinguished from the real "multi-photon" exchange which lead to inelastic transmission 
^. The power of our method is demonstrated using the example of the oscillating delta 
barrier which has been extensively studied in the literature. We show that it is possible 
to calculate and understand transmission properties which have not been accessible so far. 
An example is the transmission zero associated with the Fano resonance characterizing the 
quantum dynamics in this system. The virtual "multi-photon" processes, mentioned above, 
are relevant for the determination of its location. 
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The paper is organized as follows. In section 2 we present the main idea of the perturbative 
scheme used to describe the transmission properties of a time- dependent quantum device. 
In section 3 we present our method focusing on the case of periodic driving and using as 
a simple example the delta barrier with a harmonically oscillating strength. In section 4 
we give the results for the transmission properties of the oscillating delta-barrier setup. In 
particular we demonstrate how one can systematically calculate and interpret details of the 
behaviour of the transmission coefficient as a function of the incoming energy within our 
approach. Section 5 contains our concluding remarks. Finally, extensive formulas and their 
derivation are provided in the appendices. 



II. FIELD THEORETICAL PERTURBATIVE APPROACH FOR TRANSMIS- 
SION IN TIME-DEPENDENT POTENTIALS 

We consider the transmission of a quantum particle through a localized time-dependent 
potential. Our specific analysis is performed assuming a one-dimensional setup. The pro- 
posed approach can however be easily generalized to higher dimensional cases. Initially, i.e. 
for if: — !■ — oo, the wave function of the quantum particle for x — )■ — oo is a plane wave with 
energy Ei and momentum directed from the left to the right (incoming quantum particle). 
We also assume that the wave function for t — ?■ oo and a; — )■ oo is a plane wave with energy 
Ef and the same direction of the momentum (outgoing quantum particle). The amplitude 
for the scattering of the considered quantum particle by a specific local time-dependent 
potential is in general given in terms of the corresponding causal Green's function A„ ^ as: 

n,m 

where Pi and pf are the input and output momenta of the particle while Ei and Ef are the 
corresponding energies. In eq. ([1]) we have introduced the temporary basis {\n{t))} with n 
being a collective index for all the quantum numbers needed to fully identify each temporary 
eigenstate of the time-dependent Hamiltonian H{t) at time t. In this basis the instantaneous 
Schrodinger equation reads: 

H{t)\n{t)) = E^{t)\n{t)) (2) 
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The relevant quantity for determining the Green's function A„^m are the geometrical phases 
{n{t)\ihdt\m{t)) given as: 

; /M.^oi / ^^ {n{t)\ihdtH {t)\m{t)) / / N /„N 

{n{t)\indt\m{t)) = ^ > ('^T^'^) (3) 

The diagonal elements 7n(i) = (ri(t)|i^9t|n(t)) are the Berry phases. In the temporary basis 
the propagating kernel A„ ^ fulfils the Green's equation: 

^[{ihdt - Enit))5n,m " $n,m(t)] Am,n' = SnySit - t') (4) 
m 

where the infinite dimensional matrix $ contains only the non-diagonal geometrical phases: 

^n,m{t) = {n{t)\ihdt\m{t)) - 5n,mlm{t) (5) 

while the Berry phases are included in the effective energies En{t): 

(t) =E^{t)-J2 ^n,mlm (t) (6) 
m 

In operator notation the propagator An,m can be written as: 

An,m{tf,ti) = Htf)\-. ^= ^|m(ti)) ; Ho = -ihdt (7) 

Hq - E - ^ 

The non-diagonal matrix $ in ([7]) renders the considered problem analytically intractable 
in its general form. To proceed with the calculation of the amplitude ([1]) it is necessary to 
develop a scheme allowing for the expansion of the propagator in terms of simpler calculable 
sub-processes. The main assumption in our approach is that there exists an ordering with 
respect to the magnitude of the transition amplitudes of the different dynamical processes 
taking place in the scattering off a time-dependent potential. The latter is implied by the 
energy difference between the incoming and the outgoing state as stated in eq.(l3]). The 
amplitudes of the elastic processes dominate while inelastic processes with small energy 
transfer are more probable than those with large energy transfer. This property, if valid 
and consistently applicable, suggests that the non-diagonal matrix $ can be treated as a 
perturbation and the following expansion of A„ „^ is possible: 



In eq. ([8]) the expansion breaks down in the case of zero eigenvalues of the denominator 
Ho — E requiring a special treatment. We will come back to this point later on. Using the 
expansion ([S]) we can write the amplitude ([T]) as follows: 

oo 

n,m r=l 

where the An^m term in the above sum involves i insertions of the transition matrix $ while 
Sfi,static is the zeroth order term which does not contain $. The first term in the sum on 
the r.h.s. of eq. ([9]) is: 

/ + 00 ^ 1 
dti ~\h{tim,,iAh{ti)\-^ ~\m{U)) (10) 
-oo i^i^ Hq — h Hq — h 

and the higher order terms have a similar structure as implied by the expansion ([8]). The form 
([9]) allows the use of a particle picture for the interpretation of the dynamics in scattering off 
a time-dependent potential having at the same time a simple diagrammatic interpretation: 
the amplitude Sji is decomposed in a sum of sub-processes. Each sub-process is a sequence of 
two elementary processes: the particle propagation being in a particular eigenstate and the 
transition between two states of the temporary basis. The explicit form of Sfi depends on 
the applied potential and especially on the temporary basis. The power of the diagrammatic 
representation of Sfi is that it allows the calculation of desired properties (like the energy of a 
transmission zero or a local transmission maximum) by isolating the contributing dynamical 
processes. In general the spectrum of H{t) contains both a discrete as well as a continuous 
part. Therefore the allowed elementary processes can be classified as follows: 

• continuum- continuum (c/c) transitions 

• continuum-bound (c/b) and bound-continuum (b/c) transitions 

• bound-bound (b/b) transitions 

• propagation in a continuum state 

• propagation in the bound state 

A typical sub-process contributing to Sfi is shown in Fig. 1. The curly line indicates propa- 
gation in a continuum state while the dashed line means propagation in a bound state. The 
full black circles indicate c/c transitions while the circles containing a cross indicate a c/b 
or b/c transition. 

6 



^AAJ©- - ^ — ©w\•^AA/ 



FIG. 1. A typical sub-process involving 3 transitions (continuum - bound state, bound state - 
continuum and continuum - continuum) contributing to the transition amplitude (1) at the third 
order of the proposed perturbation scheme. 

III. THE CASE OF TIME-PERIODIC POTENTIALS: A DELTA-BARRIER WITH 
OSCILLATING STRENGTH 

The calculation of the transmission properties in a time- dependent potential significantly 
simplifies if the driving is periodic. In this section we will demonstrate how the perturbative 
scheme introduced in the previous section works in practice by performing an analysis of 
the transmission properties of a monochromatic wave passing through a periodically varying 
potential. As a concrete example we consider the transmission through a delta barrier with 
oscillating strength in one dimension. As mentioned already in the introduction, our main 
purpose is to calculate and explain features of the transmission behaviour which are not easily 
accessible by other approaches like direct integration or Floquet theory. The Schrodinger 
equation of the considered problem reads: 

z;i^^ = - — ^^-^(t)5(a:)vl/(x,t) ; g{t)=g{t + T) (11) 

with 7" = 77 being the period of the oscillating 5-potential and u the associated frequency. 
For fixed t the problem becomes static and can be easily solved. Introducing the length 
scale ^Q = and the energy scale eo = hu the static version of eq. (fTT]) can be written in 

dimensionless form as follows: 

l4pj-9rS{o]ui^;9r) = eir)ui^;gr) ; r = ujt (12) 



y 

with ^ = gr = ^a/^5'(''")) ^i'^) = and u{^; gr) is the corresponding wave function. 
For a given time instant the spectrum of the Hamiltonian consists of continuum states and 
one bound state with energy = —\g1- The later exists only during the time period for 
which g^- > 0. The associated complete and orthonormal temporary basis is: 

uti^; g^) = -^{e^''^ - ^^e^'^l^') k = V2e e > 
V27r gr + 



«6(e^7.) = x/^e-5^l«l gr>Q (13) 



as given in 18 1. 



The perturbative calculation of the scattering amplitude (eq. ([T])) requires the summation 
of all the contributing sub-processes in increasing order. According to the discussion in the 
previous section the order of a term in the perturbation series is determined by the number 
of transitions (c/ c, b/c or c/b) and can be diagrammatically presented by the series of graphs 
shown up to second order in Fig. 2. Since there is only one bound state in the temporary 
spectrum there are no possible transitions between bound states. 

V\OW=^AAA/+^Ay•V\/+ V\J•^AJ•^A/+^A®--^--e^A/+ •• 

FIG. 2. The diagrammatic expansion for amplitude (1) up to the second order in the introduced 
perturbative scheme. 

To first order only processes with a single transition, necessarily of c/c type, partici- 
pate. These may be elastic or inelastic. In the inelastic case the final energy is given as: 



^/ — \AF+^ where n takes the values n = ±1,±2,.. under the restriction that + n 
remains positive. This is a consequence of the periodic driving jj] valid for all orders of 
the proposed perturbation expansion (real "multi-photon" exchange). Before going on with 
the explicit calculation of the transmission amplitude ([T]) it is useful to present the matrix 
elements for the building blocks of the contributing sub-processes. 



Continuum- continuum transitions 

The transition amplitude from one continuum state to another has the form: 

with 9k{T) = arctan(^). The superscript (±) is omitted in $fcfc'('^) since the r.h.s. of eq. f lT^ 
is independent of these signs. As expected the diagonal term k = k' diverges and needs to 
be regularized. We adopt the usual box regularization in order to define a regularized Berry 
phase: 

97r 1 ^ k 

{k(^\T)\^h^Ak^^\r)U = y{k^^\r)\^h^Ak^^\r)) = — 7.(r) = -^9r^^, (15) 

where ^ = Jy ^ = Sr{0) with 6r{0) being the delta function regularized on a finite volume 
V. With this choice r] in eq. ( IT^ is finite fulfilling rjV cc. 1. After including also the off 
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diagonal terms we get the general, regularized expression, for ^^^'(t): 

(^{k^^\T)\^h^r\k'^^\r))-J,{r)6{k-k')) (16) 

with 5{k - k')^kk'ir) = 0. 



rioo 



Continuum-bound state transitions 

Denoting with |fo(T)) the unique bound state of the potential (existing only for r e 
(27rm, 27rm + vr) with m = 0, ±1, ±2, ...) we can express the c/b transition as: 

= ±ii(r)\.m\k'->(r)) = -2«y|^-^JL_^e«'« (17) 

Obviously the b/c transition amplitude is the complex conjugate of $6^. 
The zero-order approximation of the scattering amplitude (II]) is given as: 

5(0) = -^e^i^f-f-^^^^\A^^\rf,r,){kf\b)m) - 

POD 

/ dkAf\rf,n){{kf\k(^\rf)){k'^+\r,)\h) + (+ O -))] (18) 

In eq. (fT5]) A^"-* is the propagator in the bound state while is the propagator in a 
continuum state. The causal form of the latter reads: 



The bound-state propagator is, in the interval < r, r' < 27r, the solution of the Green's 
equation: 

[z9,-e,(r)]Af (r,r') = -5(r-rO (20) 
The causal solution of eq. (l20|) is: 

(r, t') = t9{r - r')e-'^r> dri^dn) (21) 

while the solution obeying the periodic boundary condition A^^"* (r, r') = A^""* (r+27™, r') {n = 
0,±1,±2, ..) reads: 

00 

A[°V, = ^ J2 ~ lT^m)e-'^r'^^.^<^-^-^^r^) (22) 



m=— 00 



In obtaining the S'-matrix amplitude one has to perform the limits — )■ —00, — t- 00. 
For simplicity we will assume here gij) = gosinr. Thus a consistent treatment requires 



Tf = —Ti = limAr_^oo ^ such that g{Tf) = g{Ti) = 0. In this case the zero-order contribution 
to the scattering amphtude is trivial containing no transitions (free transmission). We will 
use in the following the notation {nt, mc, lb) for the classification of the contribution of the 
various sub-processes to the S-matrix. nt means that the considered sub-process contains 
in total n transitions. From these n transitions m are of c/c type while / are of b/c or c/b 
type. Using this notation we write: 

Sf;''^''^) = S{kf-h). (23) 

The first order term, containing a single c/c transition, becomes: 

/oo 
rfr,Af;(r„n)$.,..(n)Ai°)(r„r,) (24) 
-oo 

Obviously there are no first order processes involving a single b/c or c/b transition since 
the outgoing and incoming states belong necessarily to the continuum spectrum (positive 
energy). Eq. fl2^ can be rewritten as: 

-,(li,lc,Ob) _ ^kfki{n) ^ 



^(iM.o.) ^ ^ ^^fpm - ^Jkf + 2n) + [kf -kf) (25) 



with: 



and therefore 



■inT 



2lT 



dr 







271 

Since the renormalized c/c transition obeys 



e*"" (26) 



^kjkA^f - fci) = 

the term n = is not included in the sum of eq. ( |25|) . The calculation of the amplitude A^j^., 
introduced to describe the S'-matrix contribution of the c/c transitions, is straightforward 
and leads to the expression: 

s{n) = e{n) - {-ire{~n) , q,{n) = ^{J k' + gi - kT (27) 

with Afc^.fc.(0) = 0. The various n contributing to the sum in eq. (125|) correspond to the 
usual higher Floquet modes. As it can be seen from eq. f l7r|) the amplitude A^j^ki for the 
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inelastic continuum- continuum transitions decays exponentially with \n\ for |^ < 1 (where 
A;< = min(A;j, ki)): 

in accordance with the fast convergence of the Floquet sum observed (but not explained) 
in the analysis of the dynamics of the oscillating delta-barrier in the literature. It must 
be noticed here, that, for |^ > 1, although the approximation f l25]) does not hold, the 
perturbative expansion is still valid. The reason is that, also in this case, the contribution 
of the diagrams with increasing number of transitions decreases. Since the dimensionless 
coupling qq is given by the original coupling multiplied by a factor proportional to w"^/^ our 
treatment is necessarily non-adiabatic and becomes exact either in the weak coupling or in 
the rapid oscillations (very large oscillation frequency u) limit. One important issue to be 
noticed here is that the expansion of the S-matrix in terms of the number of transitions, 
as described above, allows the decomposition of the transmission process into elementary 
sub-processes giving a consistent meaning to our approximation procedure. The emerging 
perturbative scheme can also be understood in terms of an expansion in powers of the ratio 
of the coupling qq over magnitude of the incoming wave vector k^. As can be directly 
confirmed from eqs. (fT6l) and ( fTTl) . the c/c transitions are of order ^ while c/b or b/c 

/ \3/2 

transitions are of order f ^ j . The magnitude of the various diagrammatic contributions 
to the S-matrix is then quantified by the leading power of p, a power that increases as the 
number of transitions increases. 

In the next order (two transitions) the c/b or b/c transitions are also possible. Typically, 
the corresponding sub-process is demonstrated by the fourth diagram on the right hand side 
shown in Fig. 2. The contribution of this term to the amplitude Sfi is given as: 

-oo J — oo 

(29) 

This amplitude is non-zero only if the time variables are in the region (27rz/, vr + 2'kv). 
The needed bound-state propagator fulfils the periodic boundary condition: A[°^(r2,ri) = 
A[°^(r2 + 27rm,ri) (m = 0,±1,..) and is given by eq. fl2^ where the allowed number of 
terms which must be summed up, depends on the difference T2 — ti. These facts complicate 
the calculation of the amplitude ( 129|) . We can considerably simplify things by replacing the 
time- dependent bound-state energy, in the framework of our approximation scheme, by its 
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n] 



El 

k^ 



mean value over a period: 

^b = i^ [ dre.ir) = -^g^ (30) 

This approximation is a first order estimation, being justified in the case of very fast or very 
slow (static limit) oscillations of the potential which lead to an effective, time-independent, 
bound-state energy. Higher order corrections can be obtained by expanding the bound- 
state wave function around the effective coupling value y. Following this scheme the b/c 
transition amplitude becomes in first order: 

$,,(r) = 2^kJ^ \ ^ (31) 

Introducing now the Fourier transformations: 

oo 

n=—oo 
oo 

$6A:,(r)e-2^^'^>(^)= B,kAn)e-'^^ (32) 

n=— oo 

with: 



2vr Jo 

BbkXn) = ^ 1^ dr<l,fc,^(r)e-2*^'=-(^)e*"^ (33) 
we obtain for the amplitude Bkbiji): 

BM = '^ k_]g^/2 '^M) [1 - (-1)1 = B;,{-n) (34) 

where qk{n) is given in eq. fl27|) . It is straightforward to show that also the b/c transition am- 
plitudes Bkb{n) decay exponentially with the Floquet index n when the applied perturbation 
scheme is valid, i.e.: 

(35) 

Using the expressions (!32|) we can calculate the contribution of the c/b or b/c transitions 
corresponding to the third diagram on the right hand side of Fig. 2 to the S'-matrix as 
follows: 

,(2t,0c,26) _ ^ . Bk^kM 



^(2t,oc,2.) ^ _2^^ ^ -^^0{kf)6{kf - Jkf + 2n) + {kf o -kf) (36) 
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where 

^^^'^'(")= 5: — w^wTTv — (^^^ 

no=— oo 

and ii = ei — ef,. The small parameter rj in the denominator of (1371) is introduced in order to 
ensure convergence for Tj — t- — oo and is equivalent to the demand that the continuum Green's 
functions vanish in this limit. In the absence of rj the amplitude B^^k. possesses a pole at 
ii = no- The occurrence of the pole is, from a mathematical point of view, a result of the 
perturbative expansion. In fact, the time-dependence of the potential produces an effective 
Hamiltonian containing an infinite series of geometrical phases [19|. These terms are the 
origin of a non- vanishing imaginary part rj^ that naturally appears when higher order terms, 
involving virtual transitions between the continuum and the bound state (virtual "multi- 
photon" exchange), are summed up. As a consequence the denominator of the amplitude 
(1371) never vanishes on the real axis of the energy. At the same time these higher order terms 
shift the energy Cj by 6e. This is equivalent with a change in the effective bound state energy 
of the form: 

h ^ = ~g^, - 5e (38) 

Including the higher order corrections in ( l37I) we obtain a normalized redefinition of the 
continuum-bound- continuum (c/b/c) transition: 

Bk^kA^) B^^^^{n) - gfl_^ ^,^R 

no=— oo I ^ In 

where Z„ is a normalization factor and ef , rj^ are the corrected values. In Appendix A we 
shall present the detailed calculation of the renormalized factors indicating the significance 
of the virtual "multi-photon" exchange processes. Here it suffices to note that in the limit 
^ < 1 we have Bkfbin + nojBhkii—no) ~ Q(^g\^+'^o\+\no\+i-^ which in turn means that, at least, 
Bkfki ~ 0(5'o) (fo^ n = 0, riQ = ±1 or n = riQ = —1). On the other hand the behaviour 
of the corrected factors in the denominator of (1391) depends only weakly on n and no (see 
eq. flA4p and the discussion below this equation in Appendix A): ~ 0{gl), r]^ ~ O(fi'o)- 
As a consequence, the corrected values play an essential role only when the incoming energy 
is close to an integer value. 

The expression ( l39ll is finite and well-behaved for all the values of the incoming energy. 
However when the energy of the incoming particles differs from the effective bound-state 
energy by a positive integer the probabihty amplitude to arrive at the final state passing 
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through the bound state has a sharp maximum. The impact of this maximum on the 
transmission properties will be discussed below. To complete the calculation of the S- 
matrix with two transitions we have to calculate also the contribution of the second diagram 
on the right hand side of Fig. 2 containing two c/c transitions. After some straightforward 
steps we obtain: 



r.,,(„)=r^^V-^-'-<"-'>-^"-'" (41) 



4f = -Am T.^nfj^Sik, - ^kf + 2n) (40) 

n=— oo ' 

with: 

poo ^ 

Vi^ uJn.) = I dk y 

t,-tk + l + iQ- 

L = — 00 

and Afcfc/ given by (!27|) . Due to the integration over k the regularizing imaginary part in the 
denominator of eq. f HTj) can be taken as zero. It can be further confirmed, both numerically 
and analytically (see Appendix B), that the imaginary part of the above amplitude 

oo ^ 

QTkfkAiT') = X] Tr|^fc/fc!('^ ~ l)Akik,Q) ; ki = ^kf + 21 



Z — — oo 

kj+2l>0 



which is controlled by the pole at = + / dominates over its real part 

Akfk{n - l)Akki{l) 



/•oo _ _ 

^RTk.kAn) = Pr dk\2 
JO ,_ ^ 



, ei-ek + l 

(=— oo 

that avoids the pole (Pr denotes the principal value integration). In the case p < 1, the 
imaginary part is of order 0((?q" ''^''') which, for n = 0, / = ±1, is of order 0((?q), while the 
real part is, at least, of order O^Qq). 

Already at this general stage of the analysis one can gain some insight into the expected 
behaviour of the transmission coefficient. Whenever a propagator is contained in a diagram 
contributing to the transmission amplitude we observe the appearance of denominators 
which are in fact associated with the - -"^ - factors occurring in the expansion (11 01) . These 
denominators are mainly responsible for resonant structures in the transmission coefficient. 
The most important contributions come from diagrams involving propagators close to the 
incoming and outgoing states since in this case a pole structure in the complex energy plane 
emerges. If the propagators are bracketed between transitions which do not involve the 
incoming or the outgoing state then the pole structure is integrated out leading to a smoother 
contribution to the transmission coefficient. However, the effects of these denominators is 
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also controlled by the transition amplitudes which occur in the numerator and play the role 
of residues in the final expression for the contribution of a given diagram. This interplay 
between propagation and transition (poles and residues) determines the overall behaviour 
of the transmission coefficient. 



IV. TRANSMISSION PROPERTIES 



Having calculated the amplitudes of all sub-processes involving up to two transitions it is 
straightforward to calculate the associated S'-matrix amplitude up to this order. In general 
for scattering off a potential varying periodically with time the S'-matrix amplitude is given 
as: 



Sfi = Tui0)5ikf-ki)+Rui0)5ikj+ki)+J2 T}i{n)5{kf - kj + In) + Rj,{n)b{kj + i/A;^ + 2n) 

(42) 

where Tjj(O), Tfi{n) are respectively the elastic and inelastic transmission amplitudes while 
Rii{0), Rfiiji) are the corresponding reflection amplitudes. They can be expressed in terms 
of the calculated sub-processes as follows: 

^^^^0) = 1 - JT]^^*^^^) - u^rfc,,,(o) + ... 

I n^i I 

i?.(o) = -^s^,^,^(o)-^r_,...(o) + ... 

071-7 /^Try A-TT? 

n^O: Rj,{n) = W - ^5-.,.. W - ^^^.^.^(n) + ... (43) 



The total transmission coefficient is obtained as: 



n= — oo 

Contribution of continuum- continuum versus continuum-bound state (or b/c) transitions 
The decomposition f H3i) allows us to isolate the contribution to the transmission of the 
sub-processes involving c/b or b/c transitions and to compare with the corresponding con- 
tribution of sub-processes involving exclusively c/c transitions. This is not possible in any 
other approach. In addition, as in Floquet theory, we can also here discriminate between 
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elastic and inelastic contributions to the transmission. Let us first concentrate on the domi- 
nating elastic transmission channel. According to eq. f H3|) the elastic part of the transmission 
coefficient is given as: 

I n^j I I n^i I 

2jri^ Bt,s(no)-BM-,(-no) „ 4jr'i 

= 1 - ^ L ?«(„„) - n„ + ^» - n'^''''-"" ^ - 

TIQ 

The discussion below eqs. (15^ and (jH]) indicates that the dominant contributions in 
are due to the 3fJ-B|^;,.(0) and 53rfc.fc^(0) for all values of the incoming energy except for 
ej ^ riQ — |ef I where the 95^^.(0) dominates. Thus it is suggestive to define the relative 
weight Wo{e) of contributions to the transmission coefficient involving transitions between the 
continuum and the bound state to contributions involving exclusively transitions between 
continuum states as follows: 

In Fig. 3 we plot for illustration Wq as a function of the incoming energy e, for two different 
values of the coupling qq. We see that the processes involving c/c transitions dominate for 
the complete energy regime except of a small region around the integer value = 1. 

Existence of a single transmission zero 

The dramatic increase of the amplitude of the sub-processes involving transitions between 
continuum and the bound state around the integer value 1 of the incoming energy has 
important consequences on the behaviour of the transmission coefficient with varying energy 
(see the inset in Fig. 3). To clarify this let us analyse the elastic channel in more detail. 
In eq. (jBl) the first term corresponds to free transmission, the second to sub-processes with 
two transitions involving one c/b and one b/c transition, while the third term involves two 
c/c transitions. The numerator of the second term is of order 0((?q) while the third term 
is always of the order of OIqq) [? ]. The most interesting behaviour is associated with the 
denominator of the second term. As already discussed, the imaginary part rjl^ is of order 
0{gQ) (since it does not occur for sub- processes up to ©((^q)). Therefore whenever ef = no 
the second term becomes of order one {0{g^)) and dominates. In this case the exact value 
of rjn is crucial since the second term becomes proportional to 4r- Within our approach 
it can be determined by going to higher order terms containing more than two transitions 
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FIG. 3. The relative weight wq of the contribution of sub-processes involving continuum-bound 
state transitions to those involving exclusively continuum-continuum transitions to the elastic part 
of the transmission coefficient as a function of the incoming energy. Two values of the coupling go 
are used: 0.1 (solid line) and 0.7 (dotted line). The inset displays the elastic part of the transmission 
coefficient in an energy domain around = 1 for go = 0.7. 

involving the bound state. Since the number of transitions between continuum and bound 
state is always even, the next to leading order contains necessarily four such transitions. 
The calculation of the contribution to the S'-matrix amplitude of diagrams involving four 
transitions between continuum and the bound state is straightforward but lengthy and is not 
presented here. The final result up to the order OIqq) is obtained through the summation of 
terms originating from diagrams with four and six transitions in total as illustrated in Fig. 4 
[? ]. The interested reader can find a short description of the calculation in Appendix A. 

In the limit < 1, it is straightforward to study the behaviour of the transmission 
coefficient (jHj) in the case ef = no which also depends on the value of no. Therefore let us 
first investigate the case no = 1 for effective incoming energies < ef < 1. When 1 - ef = 
Oigo) > r/^(l) both the imaginary part of the denominator as well as the normalization 
factor in the second term of eq. (j44|) can be neglected and the elastic part of the transmission 
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FIG. 4. The diagrams contributing up to the order 0((7q) to tjq and to the shift in e^. The dots 
indicate that a sum over all allowed permutations in the 6-transition diagrams must be performed. 

is well described by the formula: 



When 1 — ef = OIqq) the contribution of the real part of the second term increases and 
being negative leads to a decrease of the transmission amplitude. The effect is maximized 
when 1 — ef becomes of the order of 0((7g). Then the denominator of the second term 
is controlled, almost exclusively, by 1]^ , which is of the same order of magnitude with the 
numerator Bk.^h{l)BbkX~^) ~ ^(S'o)- Thus, in a narrow energy region around Cj ~ 1 the 
term describing transitions between continuum and bound state dominates (in the considered 
order of perturbation theory) driving the elastic transmission to zero: Tjj(O) ~ (for more 
details see Appendix A). Following a similar reasoning one can show that also the inelastic 
transmission as well as reflection amplitude become zero in the considered case while the 
elastic reflection |i?jj(0)| ~ 1. 

When the incoming energy increases beyond the position of the transmission zero the 
expression fH5|) for the elastic transmission amplitude is again valid. However, the value of 
the transmission coefficient will slightly decrease due to the fact that the imaginary part 
of TkikX^) increases with increasing energy and contributes negatively to the transmission 
amplitude. When the incoming energy is near an integer value uq > 1, the transmission 
profile considerably differs from the one around Uq = 1. This is due to the following facts. 
Firstly the coefficient Bk^b^no) Bbkii^^o) is exponentially suppressed when no increases and, 
consequently, the numerator of the second term in eq. f l44p becomes of higher and higher 
order in qq. Secondly, r]Q{nQ) is always of order 0((7q) almost independently from the definite 
value of no. Therefore, the second term is always suppressed for rio > 1 with respect to the 
other two terms. This is an issue which has not been resolved in the previous literature. 
In jl^ it is mentioned that the residues of the higher resonances considered as poles of 




(45) 
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the transmission amplitude in the complex energy plane decrease with increasing energy, 
however no specific explanation has been provided. 
Location of the transmission zero 

As already discussed, due to the higher order corrections, the transmission zero does not 
occur for = 1 but it is slightly shifted to lower values of the incoming energy. In fact, 
the exact location of the transmission zero is determined by the effective energy ef . Within 
our treatment it is straightforward to obtain not only rj^ but also the shift in the effective 
energy ef as well as the normalization factor Zq appearing in eq. 









1 2 

^ + -gl 




2Txi 




^ 1 - — 



with: 



^ 271" 

f^irio) = Bk^bino + l)Bbk,{-no - I) 

I + Zl 

4:71 

7oK) = -^srfc^fc,(o) 

c(„„) = 2 E r &M!!!L±i)Md:iZ!£) (47) 

I Jo — ^i — I 

More details of the calculation of these quantities are given in Appendix A. It is important to 
notice here that the effect of the higher order corrections leading to the shift of the effective 
incoming energy and to the appearance of the imaginary part rjl^ in the denominator of the 
second term in eq. ( l44l) starts to show up in the immediate neighbourhood of the transmis- 
sion zero, i.e. when 1 — ef = Olg^). In the literature the transmission zero is associated with 
a Fano resonance [ol] . Its exact position however has not been obtained from first principles. 



In 10| it is estimated by a continued fractions approach. Within our treatment it is in 
principle possible (see eq. fl46p and Appendix A) to obtain the position of the transmis- 
sion zero to arbitrary accuracy using increasingly order terms in the proposed perturbation 
scheme. In addition our treatment provides a physical picture for the origin of this shift: it is 
attributed to a hierarchical sequence of sub-processes involving repetitive transitions of the 
incoming particle from the continuum to the bound state and back. These can be identified 
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with a virtual "multi-photon" exchange of the propagating particle with the driving poten- 
tial. Sub-processes with increasing number of transitions become more and more suppressed. 



Further features of the transmission coefficient 

In addition to the above analysed transmission zero there occur pronounced local peaks in 
the transmission coefficient with increasing incoming energy. A typical structure is shown 
in Fig. 5(a) where a zoom in the region around the zero of the transmission coefficient in 
the elastic channel is performed. This structure is attributed to the competition between 
the real part of the second term in eq. (jBl) and the contribution proportional to 3rfc.fc^(0) 
contained in the third term. The later possesses a cusp at ej = 1 (see Appendix B and 
Fig. 5(b)) and contributes destructively (negative sign) in the transmission amplitude while 
the former decreases rapidly as the incoming energy increases beyond the transmission zero 
value. Depending on the position of the transmission zero with respect to the value 1 where 
the discontinuity of Q^Ffc^fc. (0) is located, this competition may lead to a pronounced effect 
(the observed peak structure) or not. For the occurrence of the peak structure the distance 
between transmission zero and = 1 should be at least of the order of Qq. It must be noticed 
that the discontinuity in 53Ffc.,fc-(0) is a consequence of integrated poles in the corresponding 
S-matrix amplitudes. 
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FIG. 5. (a) The elastic part of the transmission coefficient as a function of the incoming energy for 
go = 0.7 in the region around the transmission zero and (b) Qr^.fc. (0) as a function of the incoming 
energy in the same region. 
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V. CONCLUDING REMARKS 



In the present work we have developed a perturbative scheme for the calculation of the 
S'-matrix for a scattering problem involving a time-dependent external potential. Our ap- 
proach is inspired by the treatment of scattering in quantum field theory and contains as 
a key ingredient a temporary base description of the transition amplitudes. The order of 
each term in the perturbative expansion is determined by the number of transitions between 
different states of the temporary static spectrum. We have applied this approach to the case 
of a delta function with harmonically oscillating strength. Here the temporary spectrum 
may contain also a single bound state. An interesting feature of the developed perturbation 
method is that the transitions between a continuum state of the temporary spectrum and 
the bound state depend on a different power of the coupling of the potential than the transi- 
tions between two continuum states. Therefore the perturbative approach has to be applied 
with special care in order to avoid artefacts. The proposed scheme has several advantages 
providing a tool for a detailed study of the scattering process based on the intuitive particle 
picture instead of the usual wave picture which is implied by the Floquet method or the 
direct integration of the Schrodinger equation. These advantages are here demonstrated 
by establishing a systematic procedure for the exact calculation of the transmission zero 
associated with the appearance of a Fano resonance in the considered system. In addition, 
within the proposed scheme, it is possible to show that the transition amplitudes of the in- 
elastic channels decrease exponentially with increasing order in the Floquet expansion (real 
"multi-photon" exchange) a fact that explains the appearance of only one resonant energy 
in the transmission coefficient. Furthermore our approach reveals the significance of virtual 
"multi-photon" exchange processes which are responsible for a shght shift of the location 
of the transmission zero. Finally the developed method enables also the detailed analysis 
of local peak structures occurring in the dependence of the transmission coefficient on the 
incoming energy which cannot be attributed to resonances (poles of the transmission am- 
plitudes) and are usually not investigated in the literature. A disadvantage of the proposed 
method is that it is technically more demanding than the other existing methods mentioned 
above possessing several anomalies similar to those of quantum field perturbation theory. 
However these technical complications are balanced by the deeper insight gained for the 
considered physical problem opening the perspective of a systematic decomposition of the 
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time-dependent scattering process in fundamental sub-processes building up the skeleton of 
transport in driven quantum systems. 
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Appendix A 

In this Appendix we shall briefly present the calculations that lead to the corrected values 
appearing in eq. (l39l) of the text. To begin with, we note that the need to consider higher 
order corrections can be understood in terms of the following schematic expansion: 

^ + - (Al) 



{si + 6e) - no + i {t] + 6r]) it - no + irj (5. - no + irff {si - no + irff 

It is obvious that when e^k. n^, the corrections 5e and 5ri must be taken into account. The 
first process that contains a double pole term involves in total 4 transitions, all of b/c or 
c/b type: 



S 



,(4t,0c,46) 



fi 



Bkfb {n - n^) 



■ n2} 



Si 



ni=n3- 



-Airi — 

^ Ei + n^ + irj uo 

^ . Sfc 6 (n + no) (-no) ■«/ M 

-27rz > — —2 [-a (no) - [no)\ 

^ (e, - no + 17]) 
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ei + ni + IT] 
(A2) 
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In the last expression we wrote: 

oo 

£k ~ £i — i 



a K, e.) - 2 Pr / rffe ^^^K + O^^^H-^o) (A3) 
I i 



and 



(3 (no, Ei) = J2 ^^feife (^0 + l)Bbki {-I - no) , (^k = ^Jk^^+2l^ 



(A4) 



Comparing eq. ( lA2p with the expansion (lAip we get the corrected quantities appearing in 
the text: 

ef (no) ^ ei + \g'^ + a {rio, Si) (A5) 



<^/3K,e.) (A6) 
Some comments are in order at this point. The function a corrects the pole position: 

ei = no- ^g"^ - a (no, Si) (A7) 

o 

Its leading behaviour in the limit Qo/ki < 1 is obtained for no + / = ±1 and it is almost 
independent of no, having the form: 

a ~ —g'^ai — g^a2, ai^2 > (A8) 

The non- vanishing value of /3 (see eq. ( IA4p ) drives to finite values all the terms in the 
perturbative calculation of the S-matrix. Similarly to a, its leading behaviour, being of 
order O {g^), is obtained for the combination n-o + / = ±1 and it is almost independent of 
no. This behaviour is responsible for the zero of the elastic transmission amplitude: for 

~ 1 — — a, the amplitude ^B^.i^_ (0) dominates and being real, positive and of order 
0(1) cancels the free transmission contribution driving the total amplitude to zero. For 
£j ^ n-o > 1 the phenomenon is suppressed since (0) ~ O ^fi'o'""""^'') ■ 

For the quantitative analysis of the scattering process, in the energy region in which the 
"pole" behaviour is visible, one needs a more exact value of /3 and this calls for corrections 
coming from higher order processes. The 5-transitions amplitude that contains 4 b/c or c/b 
transitions and one c/c transition, gives a higher order correction of the pole position and 
we shall not discuss it here. On the other hand the 6-transitions amplitude with 4 b/c or 
c/b and two c/c transitions contributes to t], that is to the imaginary part of -B|^fc- (n): 
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where 



(A9) 
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Bk,b (no) 



Bkfb (n + ^0 



-r 



kfk, [n) 



(AlO) 



Combining the above results with eq. flA2p we find 



Vn^P (rio) [1 + 7n (no)] 



(All) 



For the elastic channel n = the last result is apparent in eq. fH6|) . When the incoming 
energy is at the realm of the "pole" , higher order corrections of the pole residue, coming from 
first order poles in the 3-transitions amplitude with two b/c or c/b and one c/c transitions 
and in the 4-transitions amplitude with 2 b/c or c/b and 2 c/c transitions, must be taken 
into account: 



n{3t,lc,2b) r,(4i,2c,26) 
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where 
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(A13) 



Combining the results of eqs. (IA5I) . (lAllI) and (IA13I) we find the normalised contribution 
of the b/c or c/b transitions to the scattering process as given in eq. fl5^ : 



(n) 



E 



(n + no) 5, 



(A14) 



ef (no) - r^o + ^'7^ ('^o) 

The above formulas capture the non-trivial behaviour of the scattering amplitude in the 
vicinity of a "pole". For example, the zero of the elastic transmission amplitude in the 
region \l — ef\ <^ rj^ is easily verified: 
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It is not difficult to show that in the same energy region and at the same order the inelastic 
amplitudes also vanish: 



Tfi (n) = Rf, (n) 



ki Bkfb {n + l] 



1 + ^^Tk^k. (0) 



1^/1 B,^,{1) 

The term n = 1 does not contribute since -8(2) = 0. For n > 1 the coefficients 

1^/1 



R{n) = T [n) 
are negligible in comparison to -R(O) ~ 1. 
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Appendix B 

In Appendix A (and in the main text) it is extensively used that the imaginary part of 
the amplitude 



r,,,^ (0) = / dk 



oo 
/=— oo 



Ak,k (-0 Akk, (0 

6i - 6k + I + il] 



(Bl) 



dominates over its real part. Besides the numerical conffimation, we present in this Appendix 
a direct proof of this statement. Using the expression (!27|) it is straightforward to obtain: 



Air ^ 
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with: 



Qk 



= {Vk^ + l-k) ■,k = k/go) (B3) 

Observing that: (|/|) ~ (fi'o/^«)''' immediately find: '^Tk^k, (0) ~ . This estima- 

tion is valid if kf = kf + 21 > k^, i.e. if / > 0. However, negative values of / are permitted 



26 



if the incoming energy is large enough. Thus, every time the incoming energy increases by 
an integer value, the possibility of an additional negative value of I occurs. For example, if 
kf/2 = 1+0 (qq) the value / = — 1 is permitted and, in fact, it gives the dominant con- 
tribution to the sum f lB2p since ki ^ O [qq) and '^STk.^. (0) ~ O {go)- This is the technical 
explanation for the cusp appearing in Fig. Oof the text. The real part of the amplitude can 
also be calculated analytically. Using the abbreviations = Vk"^ + 1 — k, = Vk'^ + 1, 
(...) we find that: 
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The first term on the r.h.s. of the last equation is: 
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The second term has the form 
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Observing that Xk ^ go, 'tfc ~ Vs'o we immediately show that Ai ^~ (jfQ, i?i ~^ (yfQ ~^ (^g ' -^2 ^~ 
and consequently that 3fJrfc.fc. (0) ~ (^q. 
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The case kf + 21 < can be treated by using the above formulas and making the change 
ki — >■ —i\ki\ in the expressions for 3?rfe.jk.(0) + iQ^k^k^{0). The final result is the same. 
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